Electron scattering from domain walls in ferromagnetic Luttinger liquids 
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We study the properties of interacting electrons in a one-dimensional conduction band coupled to 
bulk non-coUinear ferromagnetic order. The specific form of non-coUinearity we consider is that of an 
extended domain wall. The presence of ferromagnetic order breaks spin-charge separation and the 
domain wall introduces a spin-dependent scatterer active over the length of the wall A. Both forward 
and backward scattering off the domain wall can be relevant perturbations of the Luttinger liquid 
and we discuss the possible low temperature phases. Our main finding is that backward scattering, 
while determining the ultimate low temperature physics, only becomes important at temperatures 
T/J < exp(— A/A+) with J being the magnetic exchange and A+ the backward scattering length 
scale. In physical realizations, A S> A+ and the physics will be dominated by forward scattering 
which can lead to a charge conducting but spin insulating phase. In a perturbative regime at higher 
temperatures we furthermore calculate the spin and charge densities around the domain wall and 
quantitatively discuss the interaction induced changes. 
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I. INTRODUCTION 

In a one-dimensional electron system, particle-like ex- 
citations cannot survive in the presence of interactions 
leading to a breakdown of Fermi liquid theory. In- 
stead, the excitations arc of a collective, bosonic nature 
and can be described by a universal low-energy effective 
theory, the Luttinger liquid (LL)ii"— Impurities are ex- 
pected to play a particularly important role in one di- 
mension, because electrons are not able to circumvent 
them. Studies have indeed revealed that impurities in 
many cases are relevant perturbations effectively cut- 
ting the chain and therefore impeding transport at low 
temperatures.—^— More surprisingly, however, there are 
also situations where the low-tcmperature behavior in the 
presence of multiple impurities corresponds to "healing" , 
i.e., to perfect transmissioni^iS 

One of the hallmarks of the Luttinger liquid is spin- 
charge separation. This means that the normal modes of 
the Luttinger model have either spin or charge character 
and are completely decoupled. Spin-charge separation, 
however, only holds in the case of spin degeneracy. In 
the presence of a magnetic field — which leads to spin split 
bands — the normal modes of the Luttinger model acquire 
a mixed character. This situation was first studied in 
Refs. [TTIIT2I using the Hubbard model in a magnetic field 
as a starting point. One obvious question to ask is how 
impurities affect the low temperature properties in such 
a system. Since spin and charge are no longer decoupled 
we might expect new low-energy fixed points which are 
not covered by the standard Kane and Fisher picturei^ 

Of experimental relevance is, in particular, the case of 
electrons in a quasi one-dimensional wire coupled to bulk 
ferromagnetic order. Domain walls in the ferromagnet 
then act as spatially extended magnetic impurities for the 
electrons in the wire. Such systems of coupled electronic 
and magnetic degrees of freedom have received consider- 



able interest?^ spurred, in particular, by possible appli- 
cations as magnetic domain-wall racetrack memoriesJ^ 
However, the focus has principally been on how the trans- 
port properties of free electrons behave in a ferromagnetic 
wire with a domain wall, and how these spin polarized 
currents set the domain wall itself into motion. An inter- 
esting question to ask then is whether electron-electron 
interactions can become important in such cases. Aside 
from some work on mean-field^^ and Hartree-Foc k^^d'' 
interactions, there is little consideration in the literature 
of how the electronic and magnetic behaviour of quasi 
one-dimensional ferromagnetic systems is modified in a 
strongly correlated system. 

In particular, the case of a ferromagnetic Luttinger 
liquid with a domain wall present has not yet been fully 
addressed. In the limit of an infinitely sharp domain 
wall, Pereira and Mirandapi^ have considered an effec- 
tive low temperature model containing only a spin-flip 
back scattering term. Based on this model, they have ar- 
gued that the domain wall scattering in the ultimate low 
temperature limit is the magnetic analog of the Kane- 
Fisher problem.— i^iii In other words at low temperatures 
either a spin-charge insulator or a Luttinger liquid is 
found. In addition to the spin-flip back scattering pro- 
cess considered in Ref. [l^, however, also a pure poten- 
tial (spin-independent) back scattering term is allowed 
by symmctryj ^^d'' By starting from a model for an ex- 
tended domain wall we will show that such a term is 
indeed present and can be important for the physics at 
very low temperatures. Our main focus, however, is what 
happens in the more physically relevant regime of longer 
domain walls and higher temperatures. The behaviour of 
the system in this case remains unknown, and it is that 
question we wish to address here. 

There are three possible temperature regimes in this 
problem and it is convenient to introduce here some no- 
tation for them. At high temperatures we have a "per- 
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turbative regime" where the domain wah scattering can 
be treated as a small perturbation. As we consider lower 
temperatures the perturbative treatment will break down 
due to the presence of relevant operators. At first we may 
still consider the domain wall as an extended region and 
we refer to this as the "extended regime" . This regime is 
the focus of our study. At even lower temperatures in the 
rcnormalization group flow the domain wall will become 
effectively delta function like and we can treat the rele- 
vant operators as boundary terms. This regime we refer 
to as the "sharp regime", and is the regime which has 
been previously discussed in the literature J^"— In the 
sharp regime, the spin-flip and the potential back scat- 
tering terms with scattering length A+ are the possible 
relevant perturbations determining the low-temperature 
physics. We will show, however, that if we start from 
a system with an extended domain wall with length 
A ^ A-I-, then temperatures T/J< exp(— A/A+), with J 
being the magnetic exchange, are required to enter this 
regime. This regime therefore is only accessible if one 
already starts with a very sharp domain wall, a situation 
which can possibly be realized by nanoconstrictionsJ^ 

Experimentally the construction of ferromagnetic 
chains of single magnetic atoms is already possible;^ 
In such systems the ferromagnetic order is observed to 
extend over small distances separated by regions of non- 
coUinearityj^"— It is important to note that these chains 
are assembled on some substrate and therefore cannot 
be considered as fully isolated one-dimensional systems. 
Therefore the Mermin- Wagner theorem^ which forbids 
long-range ferromagnetic order at finite temperatures in 
a strictly one-dimensional system with sufficiently short- 
range interactions, does not apply. Furthermore, the sub- 
strate has consequences for the effective spin exchange 
between the atoms. The spin exchange tensor can quite 
generally be decomposed into a symmetric and an anti- 
symmetric part. In spin chains which are part of a regular 
three-dimensional crystal the antisymmetric part is often 
forbidden by inversion symmetry. For chains of single 
magnetic atoms on a substrate, on the other hand, both 
terms are expected to be present so that non-coUinear 
spin order is a generic property of such systems. The 
presence of the substrate might also disguise the Lut- 
tinger liquid properties of the atomic wire and our model 
system might therefore be too simplistic to directly apply 
to this situation. Nevertheless, it might serve as a start- 
ing point for the investigation of more realistic models. 

Other possible candidates our model might apply to in- 
clude dilute magnetic semiconductors?^ and low temper- 
ature ferromagnetic metalsi ^^'^^ Systems where the mag- 
netic and electronic degrees of freedom belong to different 
layers would also be a possible realization. Furthermore, 
we want to point out that our analysis is also valid for 
a quantum wire with a non-uniform external magnetic 
field applied. 

Our paper is organized as follows: In Sec. |TI] we in- 
troduce the model and derive the low-energy effective 
theory by linearizing the excitation spectrum followed by 



bosonization. In Sec. lIIII we consider the first order rcnor- 
malization group (RG) equations for the various scatter- 
ing processes and discuss the fixed points of the RG flow 
in the extended and sharp regimes. In Sec. IIVI we study 
an experimentally accessible temperature regime where 
the relevant scattering terms can be treated perturba- 
tivcly. We discuss different cases depending on the hier- 
archy of the different length scales present in the prob- 
lem and calculate the spin and charge densities around 
the domain wall. In the final section we present a brief 
summary of our results and some additional conclusions. 



II. THE MODEL 

We will consider an "s-d" like model in which the bulk 
magnetization and the conduction electrons are treated 
separately (though of course still coupled). We assume 
two timescales in the problem, a fast electronic one and 
a slow magnetic one. This allows us to answer the ques- 
tion of how the presence of the domain wall affects the 
Luttingcr liquid, forgetting the effect the motion of the 
domain wall will have on the conduction electrons. For 
the already mentioned dilute magnetic semiconductors 
and, in particular, ferromagnetic metals this model with 
separate electronic and magnetic degrees of freedom, ac- 
tive on different time scales, is a realistic starting point. 

The direction of the bulk magnetization of the wire 
can be described by a unit vector n{z) = cos[8(2;)]z -|- 
sin[0(z)]y and we consider, in particular, the case 
cos[0(z)] ~ — tanh[z/A] describing the spatial profile of a 
domain wall of length A situated at 2 = 0, this is plotted 
schematically in Fig. [TJ The magnetization is coupled to 
the conduction band electrons with a strength given by 
the exchange coupling J. We consider a screened, and 
hence short range, interaction V{z — z'). To simplify the 
presentation we consider here a spin- independent, SU{2) 
symmetric interaction. We want to point out, however, 
that the low-energy theory, Eqs. ^ and ([4]), obtained af- 
ter linearizing the excitation spectrum remains valid for 
a spin-dependent interaction Vo-o-' (z — z') as long as the 
interaction is spin conserving. 

We start from the following standard "s-d" 
Hamiltonianj2£ H = Hq + Hm + Hj: 



Ho 
Hi 



Az ^1(2 



2m 



(1) 



dz %pl{z)ii;„, {z)S„„, .n{z), 



^ I dzdz'i,liz)i^liz')V{z - z')i^„,iz')^„{z). 



ipl-{z) is the creation operator for an electron of spin a 
at a position z, fi is the chemical potential, m is the elec- 
tron mass, and summation over the spin indices {a, a'} is 
implied. Due to the incommensurate nature of the spin 
split Fermi points we gain no advantage from explicitly 
considering half or full filling and the filling factor is left 
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FIG. 1: A schematic view of the magnetization orientation in 
the wire. A is the length of the domain wall. 



general. We do exclude, however, the case of very small 
filling where the Fermi energy e p would become so small 
that a linearization of the spectrum would only be appro- 
priate at very low temperatures. Furthermore, we only 
want to consider the case where the filling in both bands 
is non-zero, i.e. we are not interested in the fully spin po- 
larized "half-metallic" case^ which would bring us back 
to an effective spinless fermion model. In the following 
we set h ~ 1 and ks ~ 1- 

In order to be able to linearize the system our first step 
must be to remove the spatially dependent, and in prin- 
ciple perhaps very large, magnetization. This is achiev- 
able by rotating the spin direction around the x-axis 
to a coUincar ferromagnetic alignment via the following 
gauge transformationi^i*^ H = U'''(2;)HU(z), Vo-(^) = 



Ul,,{z)i;,,{z), and U(z) 
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The interaction is 



left unaffected as it is SU (2) invariant, the magnetization 
is locally rotated to a Zeeman term Hj^j — > Hm, and a 
gauge potential is introduced: + Hq. Thus 

H = Ho+Hg + Hi + Hm arid 



Hg 



-^E / dzi:l{z)[Q'{z)rS,,,i:,,{z) 
8m / 

--E/d^ 

Am ^ J 



4'l{zW{z)a%,,d,^,,{z) 



H 



M 



(2) 



with 9'(z) = [A cosh(z/A)]~^. The gauge potential has 
been written in a manifestly Hermitian form. The first 
term of Hq is a pure potential scatterer, the next two 
terms describe spin-flip scattering. Without Hq this is 
simply a spin split band model,— see Fig. [2] In our model 
the gauge potential introduces extended scattering terms, 
active over the length of the domain wall, which have to 
be included and arc important for the low energy physics. 
The amplitudes of the spin-flip scattering and the poten- 
tial term are proportional to X/Xp and (A/Af)^, respec- 
tively, where Af is the Fermi wave length. We are here 
mainly interested in the case A ^ Xp. Except for very 
low temperatures the potential scattering term can then 
be safely neglected. On the other hand, in the limit of a 




FIG. 2: A schematic view of the spin split bands. One of 
the gi± processes is shown, which involve spin-flip scattering 
between left and right moving electrons. 



sharp domain wall A 0, then 9'(z) — > tt5{z), the scat- 
tering terms become boundary operators, and both the 
spin-flip and the potential scatterer have to be taken into 
account. We will discuss these issues in greater depth in 
section IIIII 

The Hamiltonian ([2]) is now amenable to the usual 
bosonization procedure;^ the flrst step of which is lin- 
earization via the ansatz ipcr{z) = e^'^'^'^tlja+iz) + 
e-^^f''^ij:^_[z), where kp„ = v^2m(// ± J/2). The + and 
— indices denote the right- and left-moving electrons, re- 
spectively. Note that if the Zeeman term is large we 
must linearize around the spin split Fermi points, lead- 
ing to the breakdown of spin-charge separation. After 
linearization we have Hq -\- Hm 
linearized Hamiltonians 



iJot + Hq\, with the 



H 



Off 



vp^l Az[i)\_{z)idz'4>a-{z) - %l^\^{z)id^'il)a+{z)]{i) 



where mvpa- = kpa- 

The interaction Hi can be decomposed into spin par- 
allel and spin perpendicular components and written as 



Ho 



Ha = 



Hi = 



dz 
dz 
dz 



2 

.9411 cr 



9i^_ 

2 ParPar ~r 2 P<yrP<7r 



gi\\c 



^ I dz ^^p„,.p„_,. + ?2±p^^p-_^ ^ (4) 

(T,r— ± 
fT,r— ± 

E / ClZ - -i^ParPo 

cr,r— ± 



Here we have suppressed the spatial indices and defined 
the local density Pcr± — i^lr±ip(j±- Note that the "g- 
ology" given here refers to the already rotated model, 
not the original physical picture. The chiral electrons of 
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this linearized model, physically speaking, have a non- 
collinear spin orientation throughout the wire. Umklapp 
processes scattering two left movers into right movers 
and vice versa are always neglected here due to the non- 
commensurate nature of the Fermi wavevectors. We can 
rescale the g2\\a term to include the gi\\„ process by re- 
defining g2\\a = g2\\a " 9i\\a with H2 H2. The final 
gi±a process, schematically shown in Fig. [21 can not be 
formulated as a density-density interaction. 

Finally we have our model to be bosonized. We intro- 
duce the chiral bosonic fields (j}^r{z)^ The vertex opera- 
tor is 

= -^e-^^'-^^) (5) 
yzna 

where a is a short distance cutoff. This leads to the 



following expression for the densities: 

Par{z) = ^dz(t)ar- (6) 

Thus we can write the quadratic part of the bosonic 
Hamiltonian, iJg, which is composed from Hq, Hm, H2, 
as a matrix equation 

Ha = j dz[5,$(z)]^Ma,$(z). (7) 

= ((/)-t-+, (/)-t-„, 0^+, (/)^_). The bosonization proce- 
dure is thus sufficient to re-express all but and Hq in 
terms of a diagonalizable quadratic bosonic Hamiltonian. 
The matrix, M, is 



M = — 



1 

8^ 



fA-KVF^ + 2^411^ 

■92||t 
.94_L 
92A. 



■92||t 
A-KVp^ + 2(/4||^ 

54_L 



54± 

■iTTVFi + 2.94||; 

92\\i 



4nv 



52_L 

9i± 

92U 

'Fl- 



2.g. 



4||i/ 



M is a real symmetric matrix and as such has real eigen- 
values. We can now diagonalize M and will do so in sev- 
eral steps to fully show the comparison with more stan- 
dard expressions. In order for our final normal modes to 
have positive velocities there is a condition on this diag- 
onalization procedure, given in the appendix. However, 
for any realistic microscopic model this condition should 
always be met, which can be shown explicitly in the case 
of the Hubbard model in a magnetic field.— 

We can now make two unitary transformations. The 
first is (/)o-± = (0cr T ()a)/V^- Note that 9a is the ad- 
joint of (pa and they satisfy [(j)a{z),Ua{z')] = iS{z — z') 
where 'H.a^z) = dz9a{z). This first rotation has the effect 
of uncoupling the two adjoint fields. The second trans- 
formation is to rotate to the spin-charge representation: 
'Pc/siz) = {(ptiz) i (j)i{z)]/\/2 (and similar for the 9{z) 
fields). The effect of these two rotations can be summa- 
rized as M = U-'^U-'^MUU with [l>(z)]^ ^ [^{z)YUU 
(see equations (|A.3[) and (jA.4[) in the appendix). If 
[$(z)]"^ = {(l)s,4>c,9s,9c) then our rotated Hamiltonian 
is defined by the matrix 



M 



Vb 















\ 



(8) 



Here Ks and Kc are the spin and charge Luttinger pa- 
rameters, Vs and Vc are the spin and charge velocities, 
and Va and Vb describe the coupling between the spin and 
charge sectors. These parameters are functions of the in- 
teraction strengths and Fermi velocities, and to lowest 



order can be calculated directly, see Eq. (jA.l[) in the ap- 
pendix. At the non interacting SU (2) symmetric point 
Kg = Kc = 1. In the case of spin degeneracy we find, 
as expected, Va ~ Vb ~ Q and the spin and charge modes 
decouple. 

The non-quadratic contributions from Hq in this rep- 
resentation are 



H 



f 

Gsf 



9f 



Hhsf 



, dze'(z)sin[V27r6ls(2)] 

27rma , 

sin[z(fcFt - kpi) + \/2t:(1)s{z)], 
9b 



(9) 



sm 



H, 



Gp 



H 



Gp 



, dze'(z)sm\V2Tr9s(z)] 
ZTTma J 

z{kFt + kFi) + V2^Mz)], 

dz[e'(z)]2a,0,(z), and 

dzgba[&'{z)f 
2^^F0,(z)]. 



1 
1 



E 



87rm 

a 

cos[2zkFa 



The first term, H^^p describes a forward scattering (up- 
per index "f " ) spin- flip (lower index "sf " ) process where a 
fcrmion is exchanged between the spin up and spin down 
bands but stays on the same side of the Fermi surface. 
The second term, Hq^^ is a backward scattering spin- 
flip term where a fermion is exchanged between the bands 
and also moves from one side of the Fermi surface to the 
other. The third contribution Hq^ is a potential (lower 
index "p"), spin conserving forward scattering process. 
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The final term, Hq^, is a spin conserving backward scat- 
tering process. This is the scatterer considered by Kane 
and Fisher which gives rise to the usual insulating fixed 
point 1^ All scattering processes are only active over the 
length of the domain wall (or, more generally speaking, 
the region of non-collinear spin order) where Q'{z) ^ 0. 
The scattering coupling constants for forward and back- 
ward spin-flip scattering from the domain wall are given 
by Qf/b ~ kp-f ± kpi. We have also introduced the bare 
potential scattering values gi,^ = 1/a for convenience. 
The non-quadratic interaction term Hi± is then 



H 



1_L 



(27ra)2 



dz cos[2v27r(/)s(z)] 



(10) 



This last term also corresponds to a backward scattering 
process. It stems, however, from the interaction Hj and 
therefore involves two fermions being scattered between 
the different bands and from one side of the Fermi surface 
to the other. These contributions have their simplest 
interpretation in terms of the spin and charge modes. 
However, the spin and charge modes are not eigenmodes 
of the model, see Eq. (jH]), and we also require Eqs. (|9]) 
and (jlOp in the appropriately rotated basis. 

The diagonalization of M introduces new velocities, 
Ui, for Hg and a set of {t/ '"^j parameters. It can be 
summarized as 



rp<p rp<p 

1 ^2 

rr<t> rpcf; 

^1 ^2 

fe fe 

rp0 rp0 



and 



(11) 



with the parameters {T^ as given in Eq. (jA.6p of the 
appendix. They are all known in terms of the previ- 
ously mentioned Luttinger parameters. We have also si- 
multaneously rescaled the fields to obtain two distinct 
eigenvalues rather than four. The final Hamiltonian is 
H = H„ + Hij 



H. 



Gsf 



^Gp where 



Hg ^ ^ I dz [id.Mz)r + mz)r] 



and. with summation over i and j implied: 



H, 



Gsf " 



9f 



(12) 



dze'{z)sui[V2^Tl'0,{z)] (13) 



Gsf 



H. 



Gp 



2nma , 

sm[z{kFt -kFi) + V2TTT^(j)j{z)], 
^ / dze'iz)suiV2^T,%iz)] 
sin[-z(fcFt + kpi) + V2TTff(f)j{z)], 

I ^z[&{z)fffd,cl>,{z), and 



h: 



Gp 



8i^E/d....[e'(.)f 



cos[2zkp, + ^/2^{Tf al^Tf)^,{z)] 



for the scattering terms and finally: 



H 



1_L 



(27ra)2 



dz cos[2V27ri;%(z)] 



-2iz{kp 



- 9iU^ 



2iz{kp 



The appropriate excitations of such an SU{2) asymmet- 
ric model have no obvious physical interpretation. This 
effective bosonic field theory lays the foundation of our 
further analysis. A similar model is found by Braunecker 
et alj^i^i in a different situation, and of course by Pereira 
and Mirandapi^ but without the H^^^ and Hi± terms 
which do not play any role for very sharp domain walls 
where the length scale A is no longer present. However 
they also neglect H^^ which does play a role in the sharp 
wall limit. Indeed it is this term which, when dominant, 
leads to the Kane and Fisher insulating fixed point. This 
will become clear in Sec. IIIII where the case of a sharp 
domain wall is obtained as a specific limit in our general 
analysis. 



III. LOW ENERGY PHYSICS 

In the generic case, we have three natural length scales 
present in the problem: A+ ~ {kp-f + kp^)^^ related to 
spin-flip backward scattering, A_ ^ {kp-f — kpi)~^ re- 
lated to spin-flip forward scattering, and the domain wall 
length A. In the limit of weak magnetization A_ — > oo 
and 2A+ Xp while, on the other hand, the limit of large 
magnetization, when one spin channel becomes frozen 
out, gives A+ w A_ w Af. It is crucial for the further 
analysis to observe that the relative importance of the 
forward and backward scattering terms is now not only 
determined by their scaling dimensions but also by the 
hierarchy of the three different length scales. Further- 
more, in the RG flow we must distinguish between the 
extended and sharp regimes. In the beginning of the RG 
flow we have an extended domain wall and the scattering 
terms can be treated as bulk terms of dimension d ~ 2. 
However, when the ultraviolet momentum cutoff A dur- 
ing the RG process becomes of the order 1/A then the 
extended domain wall will begin to look effectively point 
like, i.e. of dimension d = 1. The scattering terms then 
become boundary operators. At this stage the direction 
and rate of the flow of all of the operators can change, 
leading to the flnal low-temperature fixed points. For 
this effective d = 1 flow we must take the result of the 
d = 2 flow as the "zeroth order" coupling constants. In 
physical terms, the domain wall starts to look point-like 
if the electrons are correlated spatially over lengths much 
larger than the domain wall length, i.e., if J/T ^ A/a. 

Another important point to note is that for an ex- 
tended domain wall we usually have A 3> A+, i.e., the 
backward scattering terms are strongly oscillating over 
the length of the domain wall. This leads to very small 
bare effective backward scattering couplings which can 
be estimated as follows: The derivative of the domain 
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wall profile, 0'(^;), can be Fourier transformed leading 
to O'(fc) = tt/ cosh(7rfcA/2). The non-oscillating compo- 
nent of the backward scattering amplitude is then pro- 
portional to gbQ'{k = ±1/A+). Since Q'{k) is a function 
which is sharply peaked at fc = for long domain walls, 
we have \gb&'{k = ±1/A+)| - |,96 exp(-A/A+)| < 1. If 
backward scattering is relevant, then the effective cou- 
pling constant will grow under the RG flow as 

gf ^.g,exp(-A/A+)(^^j (14) 

where Tq is an energy scale of order J. Backward scat- 
tering will only have an appreciable effect if the initially 
small bare coupling has again become of order 1 under the 
RG flow. This requires temperatures T j J < exp(— A/A+) 
which are extremely small for many realistic situations. 
We therefore expect that forward scattering — ignored in 
previous investigations of the domain wall problem — will 
play the dominant role in these cases. 

Before discussing the various regimes any further, we 
want to derive the first order RG equations for the for- 
ward and backward scattering terms. We start by writing 
a functional integral partition function^ 

(15) 

with periodic boundary conditions in imaginary time t. 
Following the standard procedure we split the fields into 
fast, (j)^ , and slow, fields. Our fast fields are defined 
for A' < \k\, \u}\/ut < A, and the slow for \uj\/ui < A', 
with A an ultraviolet cut-off. Expanding the exponent 
in terms of Hi±, Hq^^j and Hq^ and performing the av- 
eraging over the fast modes wc then re-cxponentiate the 
expression to find the appropriate scaling equations. The 
flow is parametrized in terms of I, defined as A = Age"' 
and A' = Aoe-'-*'. 

We find for the Hi±^ term to first order that 

^ = 91^. [2 - 2[{Tt f + {Tir] ] . (16) 
at V ^ " 

=71 

This term is an irrelevant perturbation (71 > 2) for any 
realistic situation we are here interested in. In the limit 
of weak magnetization we can simplify the expression to 
find 7i = 2Ks- In this limit Kg > 1 and it becomes clear 
that the term is irrelevant. 

The same analysis is performed on the domain wall 
scattering terms. For spin-fiip back scattering we find: 

^ = 9b[d~l[{ftr + {ft? + {T?f + (T|)^]]. (17) 

= 76 

As already mentioned d is the dimension of the wall, for 
spatially extended walls this is 2, whereas it is 1 in the 
limit of an infinitely sharp wall. As usual the pertur- 
bation is relevant (irrelevant) if 7^ < d (7^ > d). In 



the limit of J we can simplify this expression to 
7ft = [(JsTs)-! + A'c]/2 consistent with Refs. [ll [H. If we 
consider the Luttinger parameters for the Hubbard model 
in the zero magnetization limit we see that this term is 
always relevant for repulsive interactions both for d = 1 
and d ~ 2 but can also become irrelevant in both cases 
for attractive interactions. In the general case described 
by Eq. (jl7p the relevance or irrelevance of backward scat- 
tering depends on the set of rotation parameters T^''^ and 
no general conclusions are possible. 

Similarly the spin-flip forward scattering equation 
stands as 

^^gf.[d-i [(Ttr + (Tir + (Tfr + (td^i ] . (is) 

=1/ 

In the limit of J — > this simplifles to 7/ = [{Kg)^^ + 
Ks]/2. Since Kg > 1, in this limit it follows that 7/ > 1. 
For the case of a sharp wall {d = 1) forward scatter- 
ing is therefore always irrelevant J^i^ In the case of an 
extended wall {d = 2), on the other hand, forward scat- 
tering is relevant in this limit if Kg < 2-|--\/3. The generic 
case described by Eq. p8)) is again very complicated to 
analyze. However, at least for simple microscopic models 
and in the limit of weak interactions where the rotation 
parameters can be calculated explicitly (see Eqs. (|A.7p . 
(jA.Sp in the appendix), we flnd that forward scattering 
is always relevant for d = 2. 

The potential back scattering equations are, for the 
two spin channels, 

^ = [d - liift ± + [fi ± Tif] ] (19) 

V ' 

=76<T 

with the plus (minus) sign applying for a =^ (a =^). 
In the limit of J —> we find 76(j = [Kc + Ks]/2. In 
the isotropic limit. Kg — > 1, this term is always relevant 
for repulsive interactions. In general however it can be 
either relevant or irrelevant. Indeed, it is also possible 
that potential backward scattering is relevant for one spin 
channel and irrelevant for the other. 

The potential forward scattering term has scaling di- 
mension X = 1. It will therefore be relevant for an ex- 
tended domain wall and marginal in the limit of a sharp 
wall. 

As an example, we present in Fig. [3] the two scaling 
dimensions for forward and backward spin-flip scattering 
off an extended wall as a function of ^ = J/2eF calculated 
for a U-V model with an on-site interaction U/t — 1.10, 
and nearest neighbor interaction V/t — 0.11, where t is 
the hopping amplitude. Here the parameters of our low- 
energy effective model are determined in a lowest order 
expansion in the interaction parameters. We want to re- 
mind the reader that by rotating back to the original 
physical model one finds that "spin-fiip forward scatter- 
ing" refers to an electron which passes through the wall 
without changing its spin. 
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FIG. 3: (Color online) 7/,;, and 7f,o- for forward and backward, 
spin flip and potential, scattering from an extended domain 
wall obtained for a specific U-V model (see text for details) 
as a function of ^ = J/2£f- The regions where the operators 
are relevant or irrelevant for d = 1 are marked. The condition 
Ks = 1 has been imposed by hand at the SU{2) symmetric 
point (J = 0), see appendix. 



A. RG flow in the extended domain wall regime 

At first we will focus on the d = 2 case. For extended 
domain walls (A 3> Xf) the bare potential scattering 
terms ^ (A/Af)^ are much smaller than the spin- flip 
scattering terms ^ (X/Xp). In the temperature range 
where it is appropriate to use the RG flow with d ~ 2 
we can therefore neglect the potential scattering terms. 
They will, however, become important for the RG flow 
with d = 1 at even lower temperatures discussed in the 
next subsection. 

Though with the physical parameters we consider in 
section lTVl we find both f//_b to be relevant, by varying the 
Vc,s, Kc,s, and Va,b parameters we can find regimes where 
forward scattering remains relevant but backscattering 
becomes irrelevant. 

We can now identify several regimes. We focus here 
first on the small magnetization limit close to half-filling 
where A_ 3> A+. For the backward scattering Hamilto- 
nian, Hq^^ Eq. ([S]), this means that we take the limit 
where kp^ ~ kp^ — > kp w tt/2. We can then simplify © 
and find 



h: 



dze'(z) sin[V27r6ls(z)] 



^'^^ 2'iTma 

cos[2zfci?] sm[^/2'K 4>c{,z)] 



(20) 



We now consider the case that this term is relevant, i.e., 
gb becomes large at low temperatures. Firstly, for 'nar- 
row enough' walls, i.e., for domain walls of the order of 
the Fermi wave length Xp when the 2k p oscillations do 
not cancel out the contributions, then in order to mini- 
mize the energy the fields become locked over the length 
of the domain wall in the values {0c = V27r(m -I- |), 6*^ = 



^{n + i)} or {0c = V2^{m + \), Os = v^(n + |)} 
for integer m, n. Therefore the domain wall becomes an 
impenetrable barrier for both charge and spin excitations 
and we find a spin charge insulator jiSi^ If we consider 
longer domain walls then, as already discussed previously, 
the effective bare coupling for backward scattering will be 
close to zero. Hence, for extended domain walls, A ^ Af, 
forward scattering is always the more important because 
A+ < A_. 

Therefore there is also a regime in which only forward 
scattering is relevant: either the long domain wall case 
or the case of irrelevant backward scattering (7^ > 2). 
Considering again the limit of small magnetization and 
a system close to half-filling, the forward scattering term 
in Eq. © simplifies to 



H, 



9f 



dze'(z)sin[V27r6'^(z)] 



'-'^f 27rma , 

cos[z(fcFt - kpi)] s\i\[\/2'K(j)s{z)]. (21) 



In order to minimize the energy in this case the fields 
become locked over the length of the domain wall in 
the values {^s = -\/27r(m + \),9s = -\/27r(n + i)} or 
{0s = •\/27r(m -I- \)t9s = \/2i^{n -I- |)} for integer m, n. 
In this scenario only the spin sector is frozen out and we 
find a CISO phase where the spin mode is gapped but 
charge excitations remain gapless.— In the physical ref- 
erence frame this is a state in which the incoming spin 
current docs not scatter from the domain wall and, af- 
ter traversing the domain wall, ends in the anti-parallel 
spin configuration with respect to the bulk. This means 
that the domain wall profile can no longer be taken as 
adiabatic. 

From the above considerations we have several possi- 
bilities. Firstly we can have backward scattering as ei- 
ther relevant or irrelevant. Secondly we must consider 
the relative length scales. If A+ < A then the backward 
scattering terms are small due to averaging over their os- 
cillations. A similar case holds for the forward scattering 
with A+ — > A_ in the preceding. For the case in which 
both the forward and backward scattering length-scales 
are shorter than the domain wall length we end up in 
the completely adiabatic limit, as one would expect, and 
the system shows Luttinger liquid properties ( "adiabatic 
LL"). Again we note that at extremely low temperatures 
we have to switch to a d = 1 RG flow and backward 
scattering processes will begin to dominate and can lead 
to insulating regimes. In principle, we can also be in 
the opposite regime when both forward and backward 
scattering length-scales are longer than the domain wall 
length in which case the scaling dimensions of the forward 
and backward scattering terms alone determine what the 
low-energy fixed point is. Note that this is possible with- 
out requiring a very sharp delta function like domain wall 
profile. In general, however, the low-temperature phase 
the system ends up in in the extended regime depends 
not only on the relevance of the operators, but also on 
the hierarchy of length scales. 
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TABLE I: The different phases of the ferromagnetic Luttinger 
liquid with an extended domain wall (d = 2, relevant spin-flip 
forward scattering) depending on the scaling dimensions of 
the spin-flip backward scattering operator and the hierarchy 
of the three length scales present in the problem. 





A ^ ^± 


A+ < A < A- 


A± < A 


2 - 76 > 


Spin and charge 


CISO 


Adiabatic 




insulator 




LL 


2 - 76 < 


CISO 


CISO 


Adiabatic 








LL 



The different behaviour in the extended regime which 
could be identified from the first order RG equations are 
summarized in table ID Finally, let us also comment on 
the case of a generic magnetization and arbitrary filling. 
In this case the analysis above stays valid, the spin and 
charge modes, however, get locked into more complicated 
spin and charge density wave states over the length of the 
domain wall. 



B. Fixed points of the d — 1 sharp domain wall 
regime 

For the case of a very sharp domain wall {d = 1) 
we have shown that forward spin-flip scattering is al- 
ways irrelevant. The forward potential scattering term 
is marginal and we will ignore it as well. In this case the 
possible phases are therefore determined by the scaling 
dimensions of the two backward scattering terms alone. 
(Naturally for a sharp wall the length-scales can play 
no further role.) Formally one can find the sharp do- 
main wall limit from Eq. ([2]) by taking the limit A ^ 0, 
which requires &'{z) Tr5{z). This leads to an effec- 
tive model where all the boundary scattering terms al- 
lowed by symmetry are present. A full description of 
the phase space for this model with different relevant 
perturbations present requires the solution to the second 
order RG equations to find the separatrix between the 
different low temperature fixed points. The second order 
equation for our model is more complicated than for the 
standard sine-Gordon model. A diagonal equation in the 
(j)i 's is not recovered and to perform any further analysis 
we would have to re-diagonalize the problem and then 
renormalize the model once again, repeating these steps 
until we reached the fixed point This is perhaps not 
totally unexpected as the scattering terms we are deal- 
ing with explicitly couple the normal modes. We leave 
the more involved second order RG analysis to a future 
work and focus here on what the first order equations 
can tell us. The flow of a similar model has already been 
analyzed by Araiijo, et alJ^^ using poor man's scaling. 
In that work they consider both spin-flip and pure poten- 
tial backscattering from a sharp domain wall, treating the 



interaction only perturbatively. They find phases domi- 
nated by the spin-flip and pure potential backscattering 
processes. In contrast Pereira and Mirandapi^ consider 
only the spin-flip backscattering term and hence cannot 
recover all possible low energy phases. 

There are three possible fixed points of the RG flow 
depending on the relative relevance of the three backscat- 
tering channels, gb, gb-]-, and gb^- The system can flow to 
a spin and charge insulator, an effectively spinless Lut- 
tinger liquid, or a spin-full Luttinger liquid. Further- 
more the spin and charge insulator itself can show differ- 
ent physical behaviour in the region of the domain wall 
depending on the relative relevance of the backscatter- 
ing operators. This behaviour will be confined to some 
region around the boundary. Firstly let us discuss the 
spin and charge insulating phase. If the pure poten- 
tial backscattering for both spin channels is relevant and 
dominates, ^ba < 7fc, then the system will flow to the 
spin and charge insulating flxed point already studied 
by Kane and Fisher— with the spin part of the system re- 
maining completely unaffected. If the spin-flip back scat- 
tering term is relevant and dominates, jb < "fba, then we 
have the fixed point analyzed by Pereira and Mirandai^ 
The spin-flip back scattering will tend to equilibriate the 
number of up and down spins. This fixed point must 
therefore correspond physically to a spin and charge in- 
sulator with a region of reduced spin polarization around 
the boundary. It is also possible to have the situation 
where jba < 76 < lbs with 'jba,7b < 1- In this case 
as the system approaches the fixed point first one spin 
channel will become insulating, however, before the sec- 
ond spin channel also becomes insulating the spin-flip 
backscattering term will tend to align the spins into the 
flrst channel. Once the spins are scattered into the first 
channel they are more likely to be scattered without a 
spin flip than back into the other spin channel. There- 
fore the flnal fixed point will likely be an insulator with a 
region of increased spin polarization around the domain 
wall. Secondly if potential back scattering is relevant 
for one channel but irrelevant for the other and spin-flip 
scattering is irrelevant as well, then we are left with an 
effective spinless Luttinger liquid, i.e., we have one insu- 
lating and one conducting channel. Finally, if all back 
scattering operators are irrelevant, then the system re- 
mains a Luttinger liquid. 

These results are summarized in table |TT1 



IV. SPIN AND CHARGE DENSITY IN THE 
PERTURBATIVE REGIME 

We now return to the case of an extended domain wall, 
Xf, where the potential scattering terms can be ig- 
nored and consider a perturbative temperature regime 
where the spin flip scattering terms can be treated per- 
turbatively. This allows us to calculate the spin and 
charge densities around the domain wall. Spin and charge 
density oscillations around impurities are not only ex- 
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TABLE II: The different phases of the ferromagnetic Lut- 
tinger liquid with an effectively sharp domain wall (d = 1, 
irrelevant forward scattering) depending on the relevance or 
irrelevance of the backward scattering operators. The physics 
of the system surrounding the domain wall for the spin and 
charge insulating phase depends upon the relative relevance 
of the backward scattering operators and is discussed fully in 
the text. 



7b < 1 


7()<T < 1 


7b > 1 


and/or 76t,7H < 1 


76i76ct > 1 


76t.7H > 1 


Spin and charge 






insulator 


Effective spinless LL 


LL 



perimentally relevant, but also provide a useful theo- 
retical tool to analyze the dominant physical scattering 
processes in general low dimensional strongly correlated 



systems 



41-44 



As we are interested in the case where both 



forward and backward scattering arc relevant perturba- 
tions, a pcrturbative treatment of the scattering terms 
will break down at low enough temperatures. We indeed 
find that the perturbative corrections increase as a power 
law in inverse temperature. In order for perturbation the- 
ory to be valid, we find that the following conditions, for 
forward and backward scattering terms respectively, have 
to be held: 



and 



T 
T* 



T 

^2 



T 

^2 



AttXttiT 



gs[Ttn + T^T^] 



47rAmT 
g,[TtTf + TtTl] 



(22) 



(23) 



Here T* ~ Ui/{2Tia) is a cutoff scale and the exponents 
are given by a, = [{Tff -f {T^f]/A and A = [{ft? + 
(T|)2]/4. Note that ai + as = 7/72 and /3i /^z = 7b/2 
so that the scaling dimensions show up in Eqs. \22\ and 
in the expected way. 
We want to present our results in the physical unro- 
tated frame. The perturbative results are, however, ob- 
tained in the rotated frame so that we have to use the 
gauge transformation U(z) once more. The spin density, 
in the physical frame, is 



S{z) 



1 



'ip\z)U^aUil}{z) (24) 



where the gauge rotated Pauli matrices are given explic- 
itly by 



cos[e(z)]cr2^ - sin[e(z)]c^^ 
sin[e(^)](T^'-hcos[e(^)]o•^ 



Using these relations, the spin densities in the physical 
frame S{z) can now be constructed from the spin den- 
sities Sq{z) = ^xjj^ {z)a-%l}{z) in the rotated frame. The 
corrections to the bulk in first order in forward and back- 
ward scattering off the domain wall are given by 



(A^o(^)) 



(A^o'(z)) 



(A5o^(z)> = - 



— (cos[%/27r6's(z)] X 
air 

[cos[(fci^^ - kFi)z -t- \/2TT(j)s{z)] 
+ cos[(A;Ft + kFi)z + V2T:(t)c{z) 

— (sin[v^6'^(z)] X 

[cos[(fcFt - kpijz + \f2Ti(j)s{z)\ 
+ cos[{kp^ + kp]^)z + \/2tt4)c{z)\\), and 
1 



(26) 



(27) 



dz(t)s{z) 



(28) 



-- — al„ cos[2xkFa - 2V7r</>o-(2;)] 
zaTT ^ — ' 



(25) 



Note that in the linearization procedure, strictly speak- 
ing, we should write iia{z) = ipo<j{z) + e^'^^'^^pa+iz) + 
g-tkF^z^_^(^^y As such there are iI)\{z)SiI)q{z) terms 
missing from the above spin densities which give the bulk 
values. We can also use a description where we absorb the 
effective magnetic field by a shift in the bosonic fields in- 
stead of linearizing around the spin split Fermi points. In 
this case, however, we must take into account curvature 
terms4^ This would reintroduce the bulk spin density 
terms explicitly in Eqs. ([Ml) to p8)) . 

In the following, we want to consider two examples. 
Example (a) corresponds to a case where A-|- <C A < A_ 
so that only forward scattering contributes. In exam- 
ple (b), on the other hand, we consider the case A_|_ ^ 
A_ < A so that forward scattering is again dominant but 
is oscillating over the length of the domain wall. More 
specifically, we consider values which are appropriate for 
Fe, \f ~ 0.367 nm, and a lattice cutoff a = Xp- We 
plot results for two sets of parameters: (a) = 0.112 eV, 
A = IOAf, and T = 25K; and (b) ^ = 2.23 cV, A = IOAf, 
and T = 50K. These parameters give us the following 
scaling terms: (a) 7/ = 1.12, and 7^ = 0.923; and (b) 
7/ = 1.09, and 7^ = 0.937. In both cases the condi- 
tions ([22| and ([23]) arc fulfilled with the ratio of r.h.s 
divided by l.h.s being of the order lO"'^ for the condition 
on backward scattering, Eq. (|23p . and of the order 10~^ 
for forward scattering, Eq. <\22\ . 

Figs.m and[5]show the spin density around the domain 
wall, normalized to the average value per conduction elec- 
tron, S' = i. For the situation where the length-scale 
of the forward scattering oscillations is larger than the 
domain wall length, case (a), the spin density profile is 
significantly altered, see Figs.|4Ka),[5ja). Such a shift will 
affect how the domain wall itself behaves in the effective 
magnetic field applied by the conduction electrons and 
therefore will strongly affect the domain wall dynamics. 
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FIG. 4; (Color online) The spin density Sz{z)/S, showing the 
zeroth order (dashed line) term and the total value (blue, solid 
line). The inset shows a comparison of the the zeroth order 
(dashed line) and first order (red, solid line) terms. Zeroth 
and first order refer to an expansion in the scattering poten- 
tial. All the spin figures are normalized to the average value 
per conduction electron, S — ^. The domain wall length is 



A = W\f with (a) 
4 = 2.23 eV and T 



0.112 eV and T = 25K, and (b) 



50K. 



As expected, the backward scattering plays no role in the 
considered temperature range. 

The asymmetric distortion of the spin density clearly 
visible in Figs.|4l^a) and[5{a) is due to the addition of anti- 
symmetric and symmetric combinations of the spin densi- 
ties. In contrast to Figs.|3Jb) and[5jb), where the Friedel 
oscillations are rapid when compared to the domain wall 
length, here the Friedel oscillations from forward scat- 
tering arc on a longer length-scale than A. Hence the 
changes in the spin density they cause can be seen as an 
overall distortion in its profile. 

When both A-), and A- are smaller than the domain 
wall length, see Figs. Hl^b), ^h), oscillations within the 
overall domain wall profile are clearly visible. Here the 
long wavelength oscillations are caused by forward scat- 
tering while the much faster oscillating backward scat- 
tering term causes the small "wiggles" on top of the os- 
cillations. In this case the overall spin density profile is 
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FIG. 5: (Color online) Same as Fig.|4]but for the spin density 

Sy{z)IS. 



not shifted. 

The first order charge density correction, derived sim- 
ilarly, is given by 



(Ap(z)) = e _9,0,(z) 



(29) 



Ecos[2zfcFo- - 2V7r^CT(z)]y 
/ 



The main contribution to the charge density can be cal- 
culated analytically and can be found in the appendix, 
Eq. (jA.ll[) . Results for the cases (a) and (b) are shown 
in Fig. [6] In case (a) we do see a charge build up, 
and respectively depletion, antisymmetric with respect 
to the center of the domain wall. The small "wiggles" 
on top of the overall charge rearrangement arc caused by 
backward scattering. In case (b) we see strong oscilla- 
tions of the charge density caused by forward scattering 
which are largest at the center of the domain wall. The 
very fast oscillations, which can be seen in more detail 
in Fig. [3 originate from back scattering. Finally, we 
note that the overall charge in the model is conserved, 
i.e., J dzAp(z) — 0. The total spin (5^) in the system 
is also conserved and is related to the charge density: 
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FIG. 6: The charge density correction, Ap(z), per electron 
per unit cell: po = 2e/aCm~^. The domain wall length is 
A = IOAf with parameters for panels (a) and (b) as given in 
Fig.H 



{S^{z)) = j^{p{z)). Therefore the total spin is redis- 
tributed throughout the wire in precisely the same way 
as the charge. However none of , 5*^, or S*^ arc them- 
selves conserved. 

If we conrpare Fig. [6] to the results found by Dugaev 
et al., Ref. the most striking difference is the pres- 
ence of strong Friedel oscillations. Such oscillations are 
well known and are already present in the non-interacting 
case. However, we also see that depending on parame- 
ters the amplitude of the Friedel oscillations can be quite 
small (see e.g. Fig. IHl^a)) and might be easily overlooked 
if one mainly concentrates on the overall shape. Com- 
pared to the non-interacting case^ the spin density cor- 
rections have a different profile around the domain wall 
in the transverse direction. We want to stress that our 
model offers a method of calculating non-perturbatively 
the full effect of short range electron-electron interactions 
on any physical property we are interested in. Provided 
of course that one is in the regime of sufficiently high 
temperatures such that the perturbative analysis of the 
gauge potential remains valid. This includes, in partic- 
ular, the physically important regime for quantum wires 
in ferromagnetic materials discussed above. In such sys- 



Q. 

<l 




FIG. 7: Same as Fig. ^h) where we have zoomed in on a 
shorter length making the oscillations due to back-scattering 
clearly visible. 



tems, however, the phases of the extended domain wall 
regime (see table H]) should also be accessible. 

We now want to briefly discuss consequences of the 
strong correlations in the electronic systems for the dy- 
namics of the domain wall. The magnetization dynamics 
are described by the Landau-Lifschitz-Gilbert equation, 
or some suitable generalization thereof.—"— There are 
two different aspects to this we must consider. One is the 
straightforward point that the dynamics, over the length 
of the wall, will be affected by the different spin density 
of the Luttinger liquid compared to the Fermi liquid or 
non-interacting case. The second, more interesting point, 
is whether the derivation of the non-adiabatic terms in 
the LLG equation are valid for a Luttinger liquid. 

Following Zhang and Li^ one can derive contribu- 
tions to the magnetization dynamics which allow for 
the fact that the electrons do not instantaneously fol- 
low the magnetization profile. One first writes a con- 
tinuity equation for the spins, assuming part of it to 
be always parallel to the bulk magnetization and al- 
lowing a small deviation from this. In order to derive 
the current dependent (so called /?-) terms, those which 
drive the domain wall along the wire, one assumes that 
j^'y''{z,t) = -t,BPj^'^;%z,t)n-'^''_{z,t)/e. Uz,t) is the 
charge current and P is the magnitude of the polariza- 
tion, whilst js{z,t) is the spin current. A quite reason- 
able assumption in a Fermi liquid, this of course starts 
to look more dubious in the case of a Luttinger liquid. 
In the standard Luttinger liquid model spin and charge 
are of course uncorrelated and possess different velocities. 
Thus this assumption would completely fail. For us the 
situation is not so simple as we do not have spin-charge 
separation, nonetheless what is obvious from our model 
is that spin and charge are not fully correlated. One is 
forced to work with the spin current and not the electric 
current and, as we have already seen, the spin degrees of 
freedom can behave rather differently for this model. 
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V. CONCLUSION 

We have investigated a Luttinger liquid coupled to a 
non-coUinear ferromagnetic magnetization profile in the 
shape of a domain wall. The domain wall acts as a spa- 
tially extended magnetic impurity for the electrons and 
introduces both forward and backward scattering terms, 
active over the length of the domain wall. In contrast to 
the well studied case of point-like impurities in Luttinger 
liquids the finite extent of the domain wall introduces 
a whole new layer of complexity to the problem. In a 
rcnormalization group treatment of the scattering terms, 
one has to distinguish between an extended regime and 
a sharp regime at low temperatures where the domain 
wall effectively becomes a delta function. In the extended 
regime the scattering terms are bulk operators while they 
become boundary operators in the sharp regime. An op- 
erator relevant in the extended regime can therefore be- 
come irrelevant in the sharp regime. 

Simplest to understand is the sharp regime. Here the 
spin-flip and potential forward scattering terms are ir- 
relevant or marginal, respectively. The low-temperature 
fixed points are then determined by the spin-flip and po- 
tential back scattering terms which both can be either 
relevant or irrelevant. If both are irrelevant, then the 
fixed point is the Luttinger liquid if potential back scat- 
tering is only relevant for one spin channel with the other 
terms being irrelevant, then the fixed point is an effec- 
tive spinless Luttinger liquid. In all other cases the fixed- 
point will be a spin and charge insulator. Depending on 
the relative relevance of the backscattering operators, the 
region of the insulator in the vicinity of the domain wall 
can show a reduced, increased or unaffected polarization 
in comparison to the bulk. 

If we start with an extended domain wall, however, 
then the domain wall length A will usually be large com- 
pared to the backward scattering length A+, i.e., the 
backward scattering terms will strongly oscillate over the 
length of the domain wall. This means that the effec- 
tive bare coupling — roughly proportional to the Fourier 
mode of the domain wall potential commensurate with 
the backward scattering oscillations — will be exponen- 
tially small ~ exp(— A/A_|-). So even if backward scat- 
tering is relevant, temperatures T/J < exp(— A/A+) are 
required in order to make this scattering process impor- 
tant. The sharp regime is therefore only physically rel- 
evant if we already start with a very sharp domain wall 
(of the order of a few lattice sites) which could possibly 
be realized by a nanoconstriction. 

In the extended regime, the low temperature physics of 
the model is not only determined by the relevance or ir- 
relevance of the various scattering terms but also by the 
hierarchy of the three different length scales present in 
the problem. Analyzing the first order RG equations for 
backward and forward scattering and taking the hierar- 
chy of the three length scales into account we could iden- 
tify three phases. If either, both the scattering terms are 
irrelevant, or they are relevant but the associated length 



scales arc much smaller than the domain wall length A, we 
find an adiabatic Luttinger liquid. In this case the spins 
of the electrons follow the magnetization profile and both 
charge and spin excitations are gapless. If both scatter- 
ing terms are relevant and their respective length scales 
larger than A then the first order RG equations suggest 
that the system will become a spin-charge insulator, i.e., 
the domain wall will act as a perfectly reflecting barrier. 
This case corresponds to the normal Kane-Fisher fixed 
point. Finally, there is the case of forward scattering be- 
ing relevant and the associated length scale being larger 
than A, with backward scattering being either irrelevant 
or having an associated length scale which is smaller than 
A. In this case the charge modes are gapless and charge 
is allowed to pass through the domain wall barrier. The 
spin modes, on the other hand, arc locked to a specific 
value over the length of the domain wall and the system 
becomes insulating with respect to spin transport (CISO 
phase). Physically, this means that electrons no longer 
undergo a change of spin on passing through the domain 
wall. The latter phase has not been discussed so far in 
this context and it is this phase which we believe is most 
important in possible experimental realizations. 

We also calculated the spin and charge densities 
around the domain wall for physically reasonable param- 
eters in a regime at high enough temperatures so that 
even relevant scattering terms can be treated pertur- 
batively. Here one finds spin and charge distributions 
markedly different to the previously reported mean field 
interaction case. Both the overall profile, and the lo- 
cal distribution, of spin and charge show different be- 
haviour, including Friedel oscillations. A similar result 
to ours for the lateral component of spin is found in 
the non-interacting casc^ though the transverse compo- 
nents look qualitatively different. As an outlook, we be- 
lieve that it will be interesting to study how the dynamics 
of the domain wall is changed in this temperature range 
where correlation effects dramatically alter the spin and 
charge densities compared to the non-interacting case but 
where we are still far above the phase transition temper- 
atures to the low-tcmpcrature phases discussed above. 
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Appendix: Details of the Rotations 

The spin and charge Luttinger parameters are, to low- 
est orders: 

_ VF-t + VFi , .94||t + 9iU + 254_L 



VF^; + VFi 



54|it + am - 2.94± 



Ks = 1- 



Vb 



2 47r 
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2 47r 
In order to enforce the condition Ks — 
symmetric point (J = 0) we set 

^ ^_ 1 g2||t+g2||i 
27r VF\ + f Fi 

cancehng the g2_L term by hand. It is clear that this J 
independent term must cancel when higher order correc- 
tions are included. What we do not know in this low- 
order analysis is how the J dependence of the Luttinger 
parameters will be modified by higher order terms. 
The two rotations we use on the bosonic fields are 
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To diagonalize the matrix M, Eq. we perform two 
steps. The first is a rotation to make the Hamiltonian 
diagonal, characterized by the Ri'2 terms. The second 
is a rescaling to leave us with two distinct eigenvalues 
rather than four, this introduces the F's. Together this 
gives us 



Mz) 
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Here the rotation components are 
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and 



In order for the rotated Hamiltonian to have positive 
eigenvalues the condition VcVg/v^ > K^Ks > v'^/v^Vg 
must also be satisfied. This condition seems to be always 
fulfilled, at least if one uses the integrable Hubbard model 
as the underlying microscopic lattice modelJ^ 
The rescaling of the fields requires 



and 
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Note that as the rotation is a unitary transformation we 
have the condition [R^''^]'^ + [i^j'"^]^ ~ 1. For our conve- 
nience we finally define 
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Therefore the inverse of Eq. (jA.Sp is 
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Finally, the values of the new eigenvalues are 
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As expected this reduces directly to the spin and charge 
excitation velocities in the absence of spin asymmetry. 
In such a case the Hamiltonian is already diagonal and 
the above rotation is no longer necessary. 

The analytical result for the charge density correction 

is 
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